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We provide an algebraic study of the local composite operators A^Ay — -j-A^. and A 2 ,, with d = 4 the spacetime 
dimension. We prove that these are separately renormalizable to all orders in the Landau gauge. This corre- 
^\ ■ sponds to a renormalizable decomposition of the operator A^Ay into its trace and traceless part. We present ex- 

plicit results for the relevant renormalization group functions to three loop order, accompanied with various tests 
of these results. We then develop a formalism to determine the zero temperature effective potential for the cor- 

responding condensates, and recover the already known result for (A 2 ) 4- 0, together with (a^A v — -f-A^) = 0, 
a nontrivial check that the approach is consistent with Lorentz symmetry. The formalism is such that it is readily 
generalizable to the finite temperature case, which shall allow a future analytical study of the electric-magnetic 
symmetry of the (A 2 ) condensate, which received strong evidence from recent lattice simulations by Chernodub 
and Ilgenfritz, who related their results to 3 regions in the Yang-Mills phase diagram. 



X 



I. INTRODUCTION 



— i i The dimension 2 gluon condensate (A?) in pure Yang-Mills theory has been proposed in |Q],Q|, and it has been investigated in 
different ways since then HHSHgHlslllSimillllSH. " 

In Jit] an analytical framework for studying this condensate has been developed, based on work carried out in the Gross-Neveu 
(f-) , model [15]. Different problems had to be overcome. First of all there is the gauge invariance of this condensate. In order to 
make the operator A 2 gauge invariant, one can take the minimum of its integral over the gauge orbit. Since / d xA^A^, with 
, U £ SU(N), is positive, this minimum will always exist. In a general gauge, however, the minimum is a highly nonlocal and 
thus hard to handle expression of the gauge field. A minimum is however reached in the Landau gauge (d^A^ = 0), though, so 
that working in this gauge reduces the operator to a local expression 1 . Secondly adding a source J, coupled to A 2 , makes the 
theory nonrenormalizable at the quantum level. To solve this, a term quadratic in the source must be added, which in turn spoils 
the energy interpretation for the effective action. One way around this is to perform the Legendre inversion, but this is rather 
cumbersome, especially so with a general, spacetime dependent source. One can also use a Hubbard-Stratonovich transform, 
which introduces an auxiliary field (whose interpretation is just the condensate) and eliminates the term quadratic in the source. 
Details can be found in lO]. The result was that the Yang-Mills vacuum favors a finite value for the expectation value of A 2 . The 
precise renormalization details of the procedure proposed in (3[] were given in |4|]. 



Recently, Chernodub and Ilgenfritz 111 2H have considered the asymmetry in the dimension two condensate. They performed 
lattice simulations, computing the expectation value of the electric-magnetic asymmetry in Landau gauge, which they defined 



AT c A=(^)--l-£(g 2 A?). (1) 

At zero temperature, this quantity must of course be zero due to Lorentz invariance 2 . Necessarily it cannot diverge as divergences 
at finite T are the same as for T = 0, hence this asymmetry is in principle finite and can be computed without renormalization, 
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1 We ignore the Gribov problem here, see also Ql- 

2 We shall deliberately use the term Lorentz invariance, though we shall be working in Euclidean space throughout this paper. 
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for all temperatures. The authors of lfl2Tl found that the high-temperature behavior of the asymmetry had no surprises, following 
a power law as can be guessed from general thermodynamic arguments. For the low-temperature behavior, however, one would 
expect an exponential fall-off with the lowest glueball mass in the exponent, A ~ e~'"z ]T . Instead they found an exponential with 
a mass m significantly smaller than m g \. So far, there is no explanation for this behavior. 

The goal of this work is to construct an analytical framework to investigate the electric-magnetic asymmetry studied numerically 
in 01211 . with the aim of shedding more light on the results of that paper. The hope is that something more will be found 
concerning the light mass scale m influencing the thermal behavior of this condensate. It was also noticed that the behavior of the 
asymmetry divides the Yang-Mills phase diagram into three regions in terms of the temperature. Remarkably, these regions seem 
to coincide with those associated with the condensed, liquid and gaseous states of the magnetic monopoles, whose dynamics are 
closely related to (de)confinement, see also [16]. The condensate (A^) was already related to monopoles in (US]. 

Since the computations are rather involved, we will split this work into two papers. In this paper we present the actual formalism, 
building on yllj]. Despite the fact that the quantity defined in ([]]) does not need renormalization, it is unclear how to study 
this object directly with an effective potential approach. Any finite temperature effective potential is a generalization of the 
T = potential, however the operator defining (Q]i makes little sense at T = as it is not Lorentz invariant. One could think 
about separately studying the temporal part, (Aq), and the spatial part, (A?), but this does not solve the Lorentz symmetry 
problem at T = 0: it is unclear how to couple these operators to the action such that Lorentz symmetry is maintained, while 
simultaneously allowing for a study of the ensuing effective potential. In addition, we would also like to recover the original 
results for (A*) = (A§) + (A?), i.e. we desire a clean T — ► limit. We do solve these problems by considering the operator A^Ay. 

This operator is subsequently split into its trace and its traceless part, yielding A^ and the spacetime asymmetry, A M A V — -j-A^- 
We shall prove that these two operators are renormalizable at T = 0, and that a finite effective potential can be constructed for 
these, by introducing sources J and k^, which allow to add both operators to the action without jeopardizing Lorentz symmetry. 
We also explicitly compute the effective potential at zero temperature, and we show that the only minimum is the one found 

already in yl, with a finite expectation value for A^, but with (AfjA v — j-A^.) = — as one would expect. This is already a 
nontrivial result, meaning that there is no dynamical Lorentz violation. At the end, one can take a suitable linear combination of 
both condensates to retrieve the asymmetry defined in (Q]l, at any (non)zero temperature. A priori, it is however more natural to 

study A^ and A^Ay — %-A^ than £[). In a future paper we shall then focus on the temperature dependence of both (A^) itself and 
the asymmetry. 

Summarizing the content of the current paper, we have given in section II the renormalization analysis of the operator using the 
algebraic formalism of lfl7ll . In section III we illustrate that the Hubbard-Stratonovich transformation can be used to eliminate 
terms quadratic in the action, and we compute the quantities necessary for the calculation of the effective action. In section IV 
finally the effective action itself is computed, and its minima are searched for. In section V the conclusions are presented. Some 
more technical calculations are bundled in the appendices. 



II. ALGEBRAIC ANALYSIS OF A RENORMALIZABLE OPERATOR AND ACTION 
A. Algebraic study of the renormalizability of the local composite operator A^A^ 

1. The action 

We begin by recalling the expression of the pure Yang-Mills action in the Landau gauge 

S = Sym + Sgf (2) 
= I j d d xF^ v + J d d x(b a d t AZ + <?d tl Dfc b ) , 

where 

Df^d^-gf*^. (3) 



3 



In order to study the local composite operator (LCO) A",A", we introduce it in the action by means of a BRST doublet of external 
sources (A^ v ,t1/jv), symmetrical in the Lorentz indices, in the following way 



Sk 



where s denotes the BRST nilpotent operator acting as 



CO 



CO 



2d 



J d d x (^A^+^A^c* - ^KAffi - iriASV - | V V + ^ 



(4) 



sc 



sb a = 



-D?c b , 
\gf abc c b c< , 

0, 



sK, 



/IV 



= K, 



0. 



We have shortened the notation by setting 



— K , 
TW = *\ ■ 



(5) 



(6) 



(7) 



In eq. (0]i we have used the property r^v = T| V;U , hence 

^fivAjjS^Ay) = T[ )jV A fJ D v c — T[ /jV A fl d v c . 
As is apparent from expressions (0 and ©, the action (Sym + Sgf + Sk) is BRST invariant 

s (Sym + Sgf + Sk) = . 
We notice that we can rewrite the action Sk in terms of 

(Kfn — jfiiwK) 2 = K 2 V — ^K 2 , so we might be tempted to immediately introduce a traceless tensor source k pV = K^ N — j^^K 
coupled to AfjA v . However, as not all components of k^ v can then be considered as independent due to the constraint k m = 0, using 
the derivative w.r.t. k pV becomes rather tricky, and hence also writing down suitable Ward identities. The current parametrization 
in terms of a completely general source is thus much more useful. We draw attention to the fact that we are actually coupling 
Kfj V to the (traceless) operator 



(8) 

^b/j V K, also the vacuum term — jK^K^ + ^K 2 , since 



n - -A a A"~ -^-A a A a 



(9) 



2 ^ v 2d 

meaning that we are considering the renormalization of this particular operator 0^. As we shall soon find out, we can write 
down a sufficiently powerful set of Ward identities which shall is a renormalizable operator on its 

own. 

According to the LCO philosophy ll3[ [T3l . the dimensionless parameter CO is needed to account for the divergences present in 
the vacuum Green function (o(x)0 (y)}, which shall turn out to be proportional to the specific (traceless) combination of K 2 and 
K^K^ already written down in (0]). 

Remark. We use the following definition for the derivative w.r.t. a symmetric source A^v: 

8A„ 



5A ap 2 



(10) 



2. Ward identities 



In order to translate the BRST invariance dS) into the corresponding Slavnov-Taylor identity, we introduce two further external 
sources K^, L" coupled to the non-linear BRST variations of A"^ and c" 

Sext = / d*Jt (-Kffit* + l -gf bc L a c b A , (11) 
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with 

Therefore, the complete action 

obeys the following identities 

• The Slavnov-Taylor identity 

5(E) = J d d x. 

• The Landau gauge fixing condition 

• The antighost equation 

• The ghost Ward identity 



sK£ = sL a = . 



E — SyM + Sof + Sk + S ext ■ 



8E 8E 8E 8E a 8E 8E 



0. 



8E 8E 

8c" 



g a Z = A a cl 



0. 



with 



and 



AS=^|dV" C (^X-LV) . 



Notice that the term A" p being linear in the quantum fields A a , c a , is a classical breaking. 
• Thanks to the specific way we introduced the sources K^, r\ pV and their traces K, r|, we also have 

- °' 



Sri. 



0. 



Let us also display, for further use, the quantum numbers of all fields and sources entering the action E 





A, 


c 


c 


b 


iV 




K 


L 


dimension 


1 





2 


2 


2 


2 


3 


4 


ghost number 





1 


-1 





-1 





-1 


-2 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



3. Algebraic characterization of the most general counterterm 

In order to characterize the most general local counterterm which can be freely added to all orders of perturbation theory, we 
perturb the classical action E by adding an arbitrary integrated local polynomial £ count in the fields and external sources of 
dimension bounded by four and with zero ghost number. This £ count is however restricted due to the existence of the Ward 
identities. More precisely, it amounts to impose the following conditions on £ count : 
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• The linearized Slavnov-Taylor identity 



where 



f d "" 



SiE count = q 



8E 8 8E 8 8E 8 8E 8 
^8^ + ^8^ + ^SZ? + 8L"8?' Sc fl ' """Sr^v 



obeys 



= . 



(22) 



(23) 



(24) 



• The Landau gauge fixing condition 



8E C 



8b" 



0. 



(25) 



• The antighost equation 



8E C 



8E° 



8c' 



■a T U V 



8K" 



(26) 



• The ghost Ward identity 



£ a E c ' 



0. 



(27) 



• The additional identities 





8 




8^v 




8 




Sri^v 



= 0. 



. 



(28) 



Taking into account that (/^v,!^) form a BRST doublet, from the general results on the cohomology of Yang-Mills theories it 
turns out that the external sources (A^vjT^v) can only contribute through terms which can be expressed as pure ®£-variations. 
Henceforth, the invariant local counterterm £ count can be parametrized as Il7ll 



E count = ^_j d d xF ^ F a v + q}zA -l f 



(29) 



where ao is a free parameter and A 1 is the most general local polynomial with dimension 4 and ghost number —1, given 
by 



a — 1 / Ad ( -i^a a a , t a a , 3 — a \ a , ^4 - —a—b c 

A = / d x[aiK^A p + a2L c +030^ A^ + —gf a bcC c c 

+a 5 b"-?> + a 6 ^Affi + a 7 r\A a ^ + a 8 -Tl^v + a 9 ~r[K) 



(30) 



with a\,... ,ag still arbitrary parameters. 

From the conditions (l25l l. d26l i. ( f27l > it consequently follows that 



03 = a\ 



«4 — «5 = , 



«2 = . 



(31) 



and from (I28ll we find 



«6 = 07^/, 



08 = 



(32) 
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and hence A reduces to 



A" 1 = J d d x L (k;A^' + dfjfA'V) + a 6 (^fAffi - \r^A + a 8 (%\^ - ~ ®X\k\ \ . (33) 



Therefore, for the final form of the most general possible counterterm one obtains 3 

,85 



+a 6 J d d x (^K^Affi+^A^c"- ^KA°A<<~ ^Afi.A + a 8 J d d x{^K, v K, v - ^K 2 ) . (34) 



Finally, it remains to discuss the stability of the classical action, i.e. to check whether £ count can be reabsorbed in the classical 
action E by means of a multiplicative renormalization of the coupling constant g, the parameters co, the fields {(]) ~A,c,c,b} and 
the sources {<t> = K,KfN,T\,r[fM,L,K\, namely 

E(g,co,4),4>)+#E count =E(so, <»o,4>o, 4>o) + 0(# 2 ) , (35) 
with the bare fields, sources and parameters defined as 

go = Z g g , 

0)0 = ZtoCO . 



A" 


— 7 l l 2 A a 


K 0fi - 


ZkK", , 


r " 
c 


_ 7 l/2 a 


I" - 


Z L L\ 


Eg 




K = 


ZrK , 


K 


— 1 


TlO = 


Zx\T\ , 






K(^v = 


Zk^K^ 






rio^v = 


Zr\^^v 



(36) 

and ■& the infinitesimal perturbation parameter. Notice that for consistency, we should find that Z^ (V = Z^ = Z^. 

The parameters an, a\, ag, a-/, ag, ag turn out to be related to the renormalization of the gauge coupling constant g, of the fields 
A a , c" and of the sources K, K MV , co according to 



an 



-1/2 , . „/O0 . 

Z/ = 

Z K = 1 + -r) (a 6 - a ) , 

Z(o = 1 +-r>(2a -2a 6 -a 8 ) ■ (37) 



Concerning the other fields and the sources K" LP , it can be verified that they renormalize as 



zf 


= zl' 2 = 


z b 


= z A l , 


Z K 


_ 71/2 


Z l 


_ 71/2 

— A ' 


Zn 


= Z ^v = 


Zk^ 


= z K . 



„ „l/2„l/2 , . „ / fl , "1 

ZjfZ/ Z c ' = 1 +f> [a 6 - — + — 



(38) 



This completes the proof of the multiplicative renormalizability of the LCO 0^, in the Landau gauge: the action (O is renormal- 
izable, where the Z-factor of K MV is equal to the Z-factor of K, as required. 



3 It is formally understood that we work in dimensional regularization, with d = 4 — E. We have left the d in front of the operators instead of writing 4, as this 
is important in order to get the correct finite parts once calculating in d dimensions. 
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B. Algebraic proof of the renormalizability of the local operator in combination with the LCO A 2 



1. The action 



The current problem is that the operator A 2 cannot be studied with the action (0. Indeed, if we set K^ v — Kd^ and consequently 
r) A , v = riSuv, then the action Sk = 0. As mentioned in the Introduction, from a physical point of view, we also need A 2 as operator, 
therefore we consider the following action: 

E = Sym +Sgf +Sic + Sext +Sa , (39) 

with 

= J d d x Q/A^i; + XAdc - ItA , (40) 

In this fashion, we coupled A 2 again to the action, with a doublet (X,J): 

sX = J , 

sJ = 0. (41) 

The action ( f39b allows us to study the LCOs A 2 and A , V . Clearly, these 2 operators correspond to the decomposition of AfjA v 
into its trace and traceless components. 



2. Ward identities and the counterterm 



The action E' obeys the same Ward identities as the action E (see eq. (fT3l), only the Slavnov-Taylor identity is slightly modi- 
fied: 

( 8E ' 5E ' 5E ' 5I: ' ,a^' SE' SE' SE'\ n 

5(E } - / dx [w^ + ^ +b ^ +K ^ +K -^ +J ^)=°- (42) 

and we also have an extra identity, 



r SE' r SE' 



8b 

Therefore the counterterm obeys the following linearized Slavnov-Taylor identity 

« I E' count = , (44) 

where <B% 



, i . SE 8 SE 5 SE 5 SE 8 ia 8 v 8 8 . 
Sr= /dx| 8^8^ + 5^S^ + 8^S^ + 8Z?8^ + /? S^ + ^S^+^Sl 1 ' ^ 



and it is also restricted by 

count 



hX 

After imposing all the Ward identities, we find for the counterterm 



= 0. (46) 



.-/count vcount 



J d d x (fA c ^ + J & d xa n \j 2 , (47) 
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with E count given in ( [34-b . Notice that, due to the additional Ward identity d28T l. no mixing occurs between J and K. This is a 
powerful result. A priori, a mixing between 0^ v and A^S^v cannot be excluded. 

Absorbing the counterterm ( f4Tb back into the original action E' gives the additional Z-factors (all the others are the same as 
before) 

= 1 +i3-(2ao + 2fli — an) , 
Zj = Z g Z A 1/2 = 1 --r>(ai +a ) , 

Z x = ZjZ^zl 11 . (48) 



Summarizing our result so far, we have seen that we had to introduce 2 (independent) sources to discuss the renormalization of 
A/jA v and A^, and this by means of the action E'. At the end of the story, we have identified the two independently renormalizable 
operators O^y and A 2 . This also means that the vacuum divergences ~ J 2 and the renormalization factor Zj remains unchanged 
compared with the cases already studied in JUS El]- 



III. PRELIMINARIES 

The next step will be to calculate the effective potential. Notice that this can be done relatively close along the lines of fl flltl 
We depart from the action E' and set all the external sources equal to zero, except and J. 



A. Hubbard-Stratonovich transformation 

In this section, we shall get rid of the unwanted quadratic source dependence by the introduction of 2 Hubbard-Stratonovich 
fields. We stall by writing down the complete action E', where we have set all the external sources equal to zero, except K^ v and 

E' = \ f d d xF^+ f d d x(b a d„AZ+?'d,Dfc b 



+ J d d x Q K^Affi - Y d KA "M - f W + y/ 1 ^) + J ddx {\ JA 1K - \^ 2 ) • < 49 > 

The energy functional can be written as 

e"W = y"[dA^][dc][dc][(U>]e- E '. (50) 

We now rewrite the action as 

L'=S YU + S GF + 1 d d x Q*M^o/ov- + / d d x Q/oA^A^- i^o 2 ) , (51) 



where we have used the abbreviation 



k[iv — K/jv ^—K , (52) 



and where we have used the bare quantities. We can then rewrite the action in terms of finite fields and sources. We recall that 
in d = 4 — e dimensions, we have the following mass dimensions, 

K] = — =1--, 

2 2 
4-d e 

[8\ = = 2 ' 

[J] = [K}=2, 

[Q = [(6\=d-4 = e. (53) 
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With this in mind, the action £' can be written as, 

£' = 5fm + 5 G f + J d d x 



\z K z A k pV A l ^A v - \zcZy-*atf 



flV 



d d x 



-ZjZ A JA^°~-Z^-%J 2 



(54) 



We shall now perform two Hubbard-Stratonovich transformations by multiplying expression d50t with the following uni- 
ties, 



1 = /[do]e 2Z & 



1 = / [d^ v ]e 2Z »4 



Id d x(z+±y /2 zjz A A 2 - P -^ztfy) 2 



(55) 



where we have introduced two new fields, a and cb^v We used a specific (traceless) combination, (j)^ — j-§ with (j) = (j) KK , the 
reason wherefore shall become clear soon. Let us define the following abbreviation cp^v, 



which is traceless, just like k^. Now performing the Hubbard-Stratonovich transformations yields 



,-W(J.K f , 



[dA^[dc][dc][d*][da][d<|)^]exp 



(56) 



(57) 



where 



j d d x£(A M ,c,^v) = Sym + S gf + J d d 



2Z a Z l K a> g z 2Z a Z K (ag 8Z(oCO H H 



(58) 



Notice that by the Hubbard-Stratonovich transformation, we have removed the terms ~ Xj 1 respectively ~ cofcf; v an d ~ JA 2 



respectively ~ k IJV A tl A v . The sources J and k^ are now linearly coupled to a/g and (P/jv/g as one can observe in equation d57l >. 
Hence, the usual 1PI formalism applies, with the following identification, 



(%v) = -8(0^) , 



(59) 



which follows easily from acting with ~g and on the equivalent generating functionals e L ' and e w ( J ' K vv) and setting / = 0, 



V = 0. 



Here we can also clearly appreciate the role of the traceless combination cp^v used in d55l i: it ensures that in the final action 



the traceless combination jA A ,A v — i ^f-A K A K = A , V appears, as 



1 



1 



1 



-IIVV = ^ (^v-^Mv = ~<|V ( A^y - ^-A\ \ = i U^-^M UAv-^A 2 K ) = cp„ v L) /( v. (60) 



B. Determination of the LCO parameters £ and to 



To this point, we did not yet determine the two LCO parameters £ and co. In this section we shall do so by deriving a differential 
equation for £ and CO, in an analogous way as in Q . In order to do so, we shall adapt the notation a bit, and define 8^ and 8co as 
follows, 



Co./ 2 = ^ ( i + ^ ) g 2 = /< r /:/;; • 



coo k t 



Q.jjv 



l + —)&k'=/j- t Z (a Z^k 



CO 



//V 



(61) 
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We further define the following anomalous dimensions, 



% g 2 - P(, 2 ) , 

Pfy^zj = yj(g 2 ) 

PfyhZfy = y K (g 2 ) 



(62) 



Let us first consider the case without A ,v, i.e. we set = 0. To determine a differential equation for we take the derivative 
of the first equation of ( f6Tb w.r.t. /j, 



£iC , 5C) (^+^(50) 2y,;ri ( ; - 5u : o. 



(63) 



As we can consider £ as a function of g 2 , we can rewrite the previous equation as 



8g 2 



\ii ? 2 )^(g 2 ) = 2yj(g 2 K(g 2 )+f(g 2 ), 



(64) 



with 



/(g 2 ) = e8C-Pfe 2 )^(5Q+2 Y /U 2 )5C 



(65) 



a finite function of g 2 . This particular choice of C,(g 2 ) is the unique one which ensures a linear renormalization group equation 
for the generating functional W(J) while keeping multiplicative renormalizability for (l64b is solved with a Laurent series in 



^(g 2 )=H+^+w+o(g A 



(66) 



by keeping in mind that the p-function starts at order g 4 , while a typical anomalous dimension at order g 2 . Explicit calculations 
show that f(g 2 ) at order g°. Notice also that ( |64l implies that $(g 2 ) and Jj(g 2 ) have to be known to (n + 1) loops if C,(g 2 ) is to 
be known at n loops. We refer the reader to the literature for all details involved in the LCO procedure [3, 4. [l5ll . 

We recall here from JHQj]] that we know up to three loops, 



8£ = 



N z -l 
16% 2 



35j__139 1^ (g 2 N^ | f 665 1 ( 6629 1 f 71551 [ 231 A 1^ f g 2 N } 



2 e 2 6 e / V 16^ 2 / 



6 e 3 



36 e 2 V 432 



16 



(67) 



and 



Z = 1 351 (S 2 N\ 
3 6 e V 167r 2 y 



so that 



2765 1 449 1 



72 e 2 48 £ 



\16% 2 ) 



35 / g 2 N 



16% 2 



113365 1 41579 1 



432 e j 



576 e 2 



449 fg 2 N\ 2 f 
24 ^ 16ti 2 J \ 



94363 9 
— — -— C(3) 
864 16 sv ; 



75607 



2592 16 



-77^(3) - 



^167I 2 



( g 2 N 
\16% 2 



We can now do the same for the second equation of (loTT l. Notice that from equation d62l l follows 



(68) 



(69) 



(70) 
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and therefore, 

H^~ k vN = -jKig 2 )^ ■ (71) 
Hence, we can derive an analogous differential equation for CO 

P(s 2 )^co(g 2 ) = 2y K (8 2 Mg 2 )+h(g 2 ), (72) 
as we can also make CO as a function of g 2 , where now 

h( 8 2 ) = eSco-p(g 2 ) A (8a>) +2y K (g 2 )8(0. (73) 
At the end, the generating functional W(J,K^ V ) will obey the following renormalization group equation: 

^ + $(g 2 )JL- V (g 2 ) J d 4 xJ^j-y K (g 2 ) J d 4 xVg|-)w(/,V) = 0, (74) 

which in turn ensures a linear renormalization group equation for the effective action T^a,^). There is no explicit reference 
anymore to neither K, nor co, as these LCO parameters obey their renormalization group running by construction. 



C. Loop calculations 

In this subsection we discuss the determination of the anomalous dimension of the O pV operator to three loops in the MS 
scheme, as well as the divergence structure of its associated vacuum energy. These clearly extend the analogous results for 
the original \A a2 , which were derived in flfllll. Whilst we follow the same procedures here as ifTHl . there are some novel 
features associated with treating an operator with free Lorentz indices, especially since it is also traceless. First, we note that 
we take the same general point of view as Hill and include Nf flavours of massless quarks and for the moment also consider 
the operator in an arbitrary linear covariant gauge with parameter a. At the end we will return to Yang-Mills and the Landau 
gauge a = 0. The reason for taking the more general scenario is that ironically with several additional parameters present the 
renormalization constants we will determine have more internal consistency checks. This gives us greater confidence in the 
final explicit expressions. Moreover, as the operator has free Lorentz indices, not only is there no operator mixing with the 
operator jA^ 2 , unlike that operator there is equally no mixing with any Faddeev-Popov ghost dependent operator. Therefore, the 
renormalization of the operator as defined in eq.© is multiplicative. We will use dimensional regularization in d = 4 — 2e 
dimensions to exploit the power of the calculational machinery of the MINCER algorithm, 11811 . where e = in this article. The 
algorithm is encoded, O, in the symbolic manipulation language Form, [20], and is the only viable means to achieve the goal 
of the three loop renormalization constants and anomalous dimension in a significantly reasonable amount of time. We also note 
that having the three loop results available at this stage, aside from gaining confidence in the correctness of the two loop results, 
means that the groundwork is actually laid for the future extension of our effective potential to two loops. 

First, to renormalize 0^ we insert it into a gluon 2-point function and nullify the momentum of one external gluon leg. This is 
in order to comply with the conditions of the MINCER algorithm which determines three loop scalar massless 2-point functions 
to the finite part in dimensional regularization. Concerning the gluon external leg momentum nullification, given that each 
triple gluon vertex carries a numerator momentum, this procedure does not introduce any spurious infrared singularities which 
could plague, say, a similar procedure in a scalar field theory. By contrast, the renormalization of jA'^ 2 proceeded by inserting 
the extended BRST invariant operator in a ghost 2-point function, i2ll l22tl . This reduced the number of internal gluon lines, 
resulting in a significantly fast evaluation, since MINCER relies heavily on integration by parts. An increase in internal momenta 
slows programme speed. Whilst we could insert 0^ into a ghost 2-point function, with appropriate momentum nullification, the 
three loop MINCER algorithm could only produce two loop renormalization constants for the operator renormalization constant. 
We therefore have to choose the gluon 2-point function for the operator renormalization and accept slower run times. Either 
way, however, one would still first have to give the MINCER routine scalar integrals. For the Green function we consider, 
(^a(p) o /jv{—p)Ap(0)), there are four free Lorentz indices and therefore one first needs to construct a Lorentz tensor projector. 
Ordinarily one would derive the most general tensor objects built from the momentum p** and the metric tensor, &" v , with four 
free indices in such a way that the tensor is consistent with the tracelessness and symmetry of the free indices of the operator. 
Clearly there will be more than one such independent tensor. Given whatever choice of basis is made, the scalar amplitude 
associated with each independent tensor can be deduced by inverting the tensor basis. However, as we are not actually interested 
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in the finite parts of the Green function (A^(p)o flV (—p)Ap(0)) but only its divergence structure we need only find one projector. 
This should be chosen to give a non-trivial tree contribution to ensure there is a counterterm available to absorb the infinity which 
defines the operator renormalization. Moreover, it ought to be chosen so as not to increase computer run times, especially for 
the three loop Feynman diagrams. Given these considerations we have projected with the tensor 



\uv|op 



(P) 



§up§vo 



25. 



OP" 



(75) 



The Feynman diagrams contributing to {A%{p)o )jN (— p)Ap(O)) are generated by the Qgraf package, [23], and the output con- 
verted to Form input notation by appending colour, Lorentz and spinor indices as well as the internal momenta consistent with 
the MINCER topology definitions, [20]. For this calculation there are 2 one loop, 42 two loop and 1023 three loop Feynman 
diagrams to be calculated. We follow the algorithm of [24] to extract the renormalization constant associated with A , V defined 
by 



(76) 



where o denotes a bare object. This algorithm, 02411 . in essence is such that one computes the Green functions for bare parameters 
and then scales to renormalized quantities at the end via the usual definitions. The overall divergence which remains after the 
external fields have been renormalized is absorbed by the as yet undetermined renormalization constant Z° . As this renormal- 
ization is multiplicative then this is easy to implement automatically within a Form programme. The result of our computation 
is the MS expression for the anomalous dimension 



29 g 2 N 389 fg 2 N\ 2754^(3) - 196111 / g 2 N Y 
l2 167? + 48~ VT6712J 5184 v. 16k 2 J 



(77) 



with corresponding renormalization factor 



291 fg 2 N\ 
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I 167I 2 J 
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432 e 3 



864 e 2 



7776 
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However, we have given the full expression for arbitrary gauge parameter and quark flavours in Appendix [A] For that result 
and by implication for d77l ). the anomalous dimension passes all the usual internal consistency checks. Specifically the extracted 
MS renormalization constant is a Laurent series in e. Therefore the double and triple poles in £ are predicted by the simple 
poles from lower loop order due to the renormalization group equation since we are dealing with a renormalizable operator. Our 
renormalization constant successfully satisfied this check. As an additional check on our computer code, such as the correct 
Qgraf generation of diagrams and mapping to Mincer and Form syntax, we have rerun the complete code for the calculation 
again but instead used the operator \A^ 2 . We obtained the known three loop anomalous dimension for the Landau gauge, 



lull 12211 . The arbitrary a expression in this instance is not actually a check for this purpose for several reasons. First, the 
arbitrary a result which is available is for the extension of the operator to the BRST invariant object and only for the non-linear 
arbitrary covariant gauge known as the Curci-Ferrari gauge, 12511 . rather than a linear covariant gauge. Aside from this, if one 
merely considered jA^ 2 on its own in an abritrary linear covariant gauge then there is a potential mixing to the Faddeev-Popov 
ghost mass operator c a c a for a + 0. It is not necessary to investigate this as we are only primarily interested in checking a 
working code whose output satisfies several stringent checks already. 



The second main three loop result which we have determined is the divergence structure of the associated operator vacuum 
energy. As outlined in Hill this can be deduced by considering a massless 2-point function with the operator present at each 
of the two external points and a non-zero momentum flowing in one operator and out the other. The divergence structure of 
this Green function reproduces the divergence which occurs in the vacuum graphs composing the effective potential. In essence 
these vacuum diagrams involve the constant current of the LCO formalism which gives massive propagators. By formally 
differentiating with respect to these masses one accesses the part of the vacuum diagrams which are divergent. This divergence 
can then be extracted by formally setting the mass to zero in these cut open vacuum graphs to produce the massless 2-point 
function we compute, 111 111 . As the set-up we have described involves massless 2-point Feynman diagrams, it is naturally 
accessible to the MINCER algorithm. Moreover, one does not need to nullify any external momenta. Therefore, we have 
generated the diagrams using Qgraf, which gives 1 one loop, 7 two loop and 127 three loop Feynman diagrams. However, 
as with the renormalization of 0^ we need to project with an appropriate Lorentz tensor to have scalar amplitudes, since our 



Green function is in effect 
projector 



(Of lV (p)o ap (— p)). Again only the divergence structure is required and therefore we applied the 
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(79) 
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Consequently we find the divergence is 
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(80) 



Again this is for Yang-Mills and we have devolved the full Nf expression for the Landau gauge to AppendixlAl Moreover, there 
we also give the arbitrary Nf results for various quantities derived from these expressions whose Yang-Mills versions appear 
subsequently in the main text hereafter. 



D. Solving the renormalization group equations for i^(g 2 ) and a(g 2 



For these calculations, we shall require the two loop p-function, 

P(g 2 ) = -^ 2 -2(pog 4 + Pig 6 + 0(y 8 )) , 



with 
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(82) 



In fact, in fl fTlll . S£, and were already calculated, but for the benefit of the reader, we shall repeat the results here in a 
structured way, especially since we have used a slightly different notation than in 01. 

Firstly, combining jj and 8^ into expression (l65l l yields, 



f(g 2 



(N 2 -l) 

167C 2 



(83) 



and analogously, we find h(g 2 ) by combining Jk and 8(0 into expression d73l . 
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(84) 



Secondly, solving the differential equation d64l > and invoking expression d69l and (l83l we can determine X, to one loop or- 
der 4 , 



s = 



N 2 -l 

167C 2 



9 167C 2 161 



13 g 2 JV 52 



Likewise, solving the differential equation d72l and using d77l ) and d84l i, we can determine co to one loop order, 



co = 



N L -\ 



167C 2 



1 167C 2 



73 



4 g 2 N 1044 



Thirdly, from expression ( |6TT >, we can determine Z 2 Z^ 



and ZgZ a , 



1 + 5 » = 1 _ 1 /^ 

co 3e V 16ji 2 



(85) 



(86) 



(87) 



(88) 



4 In principle, we can go one loop further with the results known. However, we shall not need this next order, as we shall only determine the effective potential 
to one loop order. It will become clear that this is already a highly complicated task. 
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IV. THE EFFECTIVE POTENTIAL AT ZERO TEMPERATURE 



A. Computation of the effective potential 



We are now ready to calculate the one loop effective potential. Firstly, from the results in section lHIDI we find that the tree level 
mass associated to A^A^ is given by 



m 



13 N 
~9~N 2 -l 



eg = o'g , 



where we have defined 



0=0- 



13 N 



9 N 2 - 1 ' 

Analogously, the tree level mass matrix associated with A^Ay is given by, 



M„ v = 4 



N 



N 2 -l 



where we have defined 



N 



N 2 -l 



(89) 



(90) 



(91) 



(92) 



To calculate the effective potential for a', <p', we can rely on the background formalism. At one loop, only the integration over 
the gluon field gives a nontrivial contribution, and we find 



1 -r(')(a',( P ' ) _ 
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where trln Q^ v is defined as, 



trln Q^ v = trln ( -3 2 8 jUV + ( 1 - j ) + S^m 2 +M^ V . . 



(93) 



(94) 



with Muv the traceless matrix which describes the asymmetry in d dimensions. 

Let us first start with the calculation of d94b . We can parameterize the traceless matrix M pV as follows, 

r 
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cl-l 
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(95) 



as we can assume that there is still spatial symmetry when separating the temporal component. We can now rewrite expression 



trln 5 



trln Q^ v = trln ^-9 2 5 A , V + f 1 - |J d p d v + S^r, 
which becomes in the Landau gauge limit £, — ► 0, 

trlnQ^v = {d- l)trln(-3 2 + m 2 )+tiin(-3 2 )+trln ( S^v 



d 2 +m 2 



a 2 +m 2 v" 1 a 2 



-8 2 + ^, 



M Xv , (96) 



(97) 
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The matrix ofthe last logarithm has the following eigenvalues: 1 — _g +m i jry with multiplicity [d— 2), 1 + _ d i +m i ^1 — J7~r 
and 1 . Therefore, we obtain 

1 A 



trlnQ^v = (,i-l)trln(-a 2 +m z )+trln(-a z ) + (rf-2)trln 1 



-d 2 + m 2 d - 1 



-trln 1 



d dl 



trln(-3 2 ) + (t/-2)trln -3 



d-l 



trln -3 2 + m 2 +A 1- 



d~y- 



iii 



Firstly, in dimensional regularization, we know that 5 

I = trln(-3 2 ) = /_— -ln/t 2 = 

Secondly, part II is a readily evaluated as 

II = (d-2)trln( -d 2 + m 2 
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(98) 
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(100) 



where we have worked in the MS scheme. Notice that, for II to be real valued, we have the following constraint: 

m 2 -- > 0. 



(101) 



Finally, part III requires much more effort due to the presence of the temporal derivatives. In the appendix, we have worked out 
this calculation in great detail. In dB30t . we have ultimately found, 
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with 2F1 a hypergeometric function. As one can find in the appendix, this part is real valued if 

m 2 +A > 0, 

which is a second constraint. Taking parts II and III together yields, 
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The remaining hypergeometric function can be expressed in terms of elementary functions, since 



5 1 1 

128z-l 7V 2 



y/z (-15z 2 +4z + 3) + I (5z 3 - 3z 2 - z - l) (ln(l + y/z) - ln(l - yfz)) 



(102) 
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(105) 



5 An overall factor (VT) will always be omitted in all the calculations. 
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In our case, we have z = V and we have 3 cases. Let us first consider z > 1, or A > 0. In this case, the In ( A ^$ m i \ picks up ; 
imaginary part, leading to 



1m 



18(471)- 



[5A 2 + l2Am 2 + 9m 4 ] In 



A — 3w z 



2887t 



(5A 2 +12Am 2 + 9m 4 ) 



Then, keeping the constraints (I10U . ( 1103b in mind and using (11051 ), we also find 
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(107) 



This is a first consistency check, as the imaginary parts neatly cancel. 

Secondly, for z < 1, or A < 0, all constituent functions in the trln are real to start with. 

A third interesting test is the z — > 1 case, i.e. the A — > limit, which should give back the known result for (d — l)trln(— d 2 +m 2 ). 
All the terms in 1 /A nicely cancel, and we recover the correct value, 
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(108) 



The second part of expression d93l is, to the orders in g 2 and e we need, readily evaluated as 
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(109) 



where we have used the results of the previous subsection lHIDI in particular d87l > and 

In summary, taking all the results together, we nicely find that all the divergences cancel in the effective potential, which, after 
some algebra and simplification leads to the final result, 
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B. A check: a handmade calculation of Sro at lowest order 



Having determined the logarithmic determinant in the previous subsection, we are now also in a position to provide a handmade 
calculation of the counterterm Sco, to have a check on the automated output leading to d80l >. First, we single out the one loop 
divergence in W(J = 0,K^ V ), as defined in (l50l >. We can identify = Mu V , meaning that it is actually given by the ~ part 
of 



N 2 -l , „ 
— 5 — trlng^v 
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Using the conventions ( f50b and ( Bil l, this divergence must be canceled by the counterterm contribution 

with the parametrization d95l l. Employing (11041 ). we find 

N 2 -l , „ 7N 2 -1A 2 c . 

— — trlne /jv = --- r7 - r — + finite, (113) 



9 16ti 2 2e 



hence, at lowest order 



Id -I N 2 - ll 7 A? 2 - II 
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a value consistent with ( f80b . 

Let us also include a check using the renormalization group equation. In particular, we have for the complete one loop functional 6 



wW(J,fy) = -~J& - - 5%/ 2 - ^i-/ 2 + ^(^V 2 - l)trlne„ v , (115) 



W ( 1 ' (/, A„v ) = - - — -CO i — + — — — ^ In - _f ^ (7 A 2 ) + terms in J, higher order terms . (116) 

2 g z 3 2 4 36 16tt/ jU z 

Focusing on the A-sector, then the application of the renormalization group equation d74T i leads to 

,|^(/,V) = ^^ (117) 

after filling in the numbers for coo, po and y^.o, since the anomalous dimension of A is determined by — jK(g 2 ), which follows 
directly from the definition of A. We have thus explicitly verified at lowest order that the definition of the LCO parameters CO 
and £ gives a generating functional W(J,K/j V ) consistent with the renormalization group equation. 



C. Minimum of the potential 



We shall now determine the minimum of the potential d 1 19b . by solving the associated gap equations 

3y(U 
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. 
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(118) 



and checking which potential extrema is a minimum. 

Before doing this, we can already check whether for small A, the minimum m 2 — (2.03A„s) 2 from (HQl]] is a stable one. A 
Taylor expansion around A = yields 
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(119) 



6 We did not write down the counterterm contribution here. 
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The first part is exactly the first order potential found in 01 [Dl, and the second part is a positive correction 7 in A 2 . Moreover, 



one can also check that 



a 2 vW = d 2 v( 

dm 2 dA dAdr 



y = at 



A = 0. 



'4.12 



(120) 



while it is already known from [IH that > at J 120b . This means that ( 1120b is a stable (local) minimum, using the 

Hessian determinant. This is a nice and nontrivial verification of the results of i3i [Tlll . while establishing Lorentz invariance (as 
A = 0), a conditio sine qua non at T = 0. 

Let us now look for other minima. The gap equations dl 181 ) yield 
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(121) 



After determining these conditions, we still need to choose an appropriate value for /j. In order to compare possible other minima 
with the minimum ( 1120b . we should operate with the same scale. Therefore, we set 



while with this choice, the coupling constant becomes, 
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■ 0.19 . 
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which is sufficiently small to assure a trustworthy perturbative expansion, as we have carried out. 

If one solves the two gap equations numerically, besides the minimum ( 11201 ), one only finds the maximum 

A = 0, m 2 =0. 

In summary, the potential dl 19b has only one minimum J 120b . 



y( 1 )( OT 2 =4.12,A = 0) 



-3.23A1 



(124) 



(125) 



where we have set = 3. As A = in the minimum, Lorentz invariance is preserved as required in this T = case. This is a 
good test for the solidness of our framework, before going to the more complicated case of finite temperature. 



V. CONCLUSION 



In this paper we have introduced an analytical formalism which allows one to study the electric-magnetic asymmetry of the 
dimension two condensate (A 2 ) first considered in [12] by means of lattice simulations. 



7 We recall that we set \? = m 2 in 13. [Till based on the renormalization group. 
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The first main result is that such a formalism exists: it is possible to compute the expection value of A 2 and of its asymmetry 
simultaneously, generalizing the formalism developed in JH. We do this in the Landau gauge, where J d^x A 2 reaches its 
minimum when A A , moves along its gauge orbit, giving a gauge invariant meaning to (A 2 ,) in the Landau gauge. Adding terms 
quadratic in the sources makes the model renormalizable, and performing a Hubbard-Stratonovich transform eliminates these 
terms, introducing new fields but restoring the usual energy interpretation of the effective action in a simple fashion. 

As a second result, we have computed the parameters which appear in the formalism. This part of the calculation is based on the 
renormalization group. Once these parameters were fixed, we have determined the effective action itself. At zero temperature it 
has the same minimum as the one already found in [3] without considering the electric-magnetic asymmetry: there is a nonzero 
value for (A 2 ), and the asymmetry is zero as must be the case based on Lorentz invariance. No other minima were found. 

Now that a consistent formalism has been developed at zero temperature, we can consider the finite temperature case. l26ll 
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APPENDIX A: HIGHER LOOP CALCULATIONS IN MOST GENERAL FORM 



In this Appendix we record the more general forms of the various LCO quantities associated with the traceless operator. In 
essence they contain the dependence on Nf massless quarks. We have checked that in the Nf — > limit they return to the Yang- 
Mills values we have used and recorded in the main body of the article. Though for the first two quantities, Zk and 8co, we have 
also included the gauge parameter dependence. We have 
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10368 
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a 



(167t 2 ) 2 
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[(3oc-29)Ca + 167>N/] 



+ [(6a 2 +43a-389)C^ + 248C A r F A^+192C F 7>A^] — 

+ [(567a 3 + 486^(3)a 2 + 3258a 2 + 3024^(3)a + 15750a + 2754^(3) - 1961 1 l)C\ 
+ (268672 - 81216^(3) -5688a)Cl7>AT / + (103680^(3) + l9l52)C A C F T F N f 



- 4U52C A T^Nf - l036SCpT F N f -25344C F T^Nf] 



2592 



+ 0(a 4 ) 



(A3) 
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[a 2 + 4a + 7 



Na 
12 



+ [(656 + 64a)T F N f - (3a 3 + 129a 2 + 287a + 1345 )C A ] 



Naci 
432 



[(7152a 2 - 22496a + 165888C(3)a - 191568 + 829440^ (3 ))C A T F N f 
- ((162a 4 + 1944a 3 + 17820a 2 + 123 12a + 239922)^(3) 

+ (585a 4 + 4698a 3 + 9768a 2 + 105418a + 266787))C] 
+ (134784a- 165888£(3)a + 509760 - 663552^(3 ))C F 7>AT / 

N A a 2 



+ (20480a + 90240)A^7 F n 



31104 
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Finally, we find the Landau gauge version of co for Nf + is 

7 



N A 

CO = 



167I 2 



4[7C A -87>JV / ] 

„2 

+ (51 1C? - 2452C A T F N f + \5\2C F T F N f + 13127iJV f 2 ) 7 ^ = = 

V A A ^ ; * * * F f, [29C A -l6T F Nf}[7C A -8T F N f }576K 2 

+ ((17352846^(3) - 10661959)C^+ (75444728 - 1 17227088^(3) )ClT F N f 

+ (123538176^(3) - 83836800)^7^ + (86994432^(3) - 62449632)C 2 T F jV / CV 

- 2 1 04704C A d. T F Nf + (59943 1 68 - 993 1 1 6 16^ (3) )C A C F T F N 2 

+ (32150528 - 35389440C(3))Ca7> 3 A^ 3 - 38502407^ 

+ (28311552t;(3)-12865536)7|A^Q ? + 81285 HT^Nfcj) 

x 



Hence we can deduce that 



(A5) 



[29C A - l6T F Nf][nC A - 8T F Nf] [7C A - 87>A^]26542087t 4 

L [HC A -47>jV / ]g 2 
24;i 2 e 

f (38857C] - 34948ClT F N f - 22ll6C A C F T F N f 

+ 6880^ + 8064^) ^[2^^Ri + ^ ' (A6> 



APPENDIX B: DETAILS CONCERNING THE CALCULATION OF PART III OF THE EFFECTIVE POTENTIAL 



In this Appendix we shall calculate step by step part III of equation (1941 . which is far from trivial. 

Ill = trlnf-a 2 + m 2 +Afl- 



d-ld 2 

- — -r In k +m +A 1 £ 

(2n) d V V ^-lA: 2 

d d fc 

(2nj d '"\~ ""V" d-l 

d d k 



r d d k 
J J2W 



In ( k A + k 2 m 2 +A\k 2 - - — T kl ) ] - trln/t 2 



^ )d ln [k^ + k^m' +2kf ~ — j +kf(ki + m 2 +A) j , (Bl) 

where the notation kj refers to (d — l)-dimensional spatial part of k. In this case, notice that for III to be real valued, we must 
have 

m 2 +A > 0, (B2) 
where we have assumed that expression (1101b is certainly fulfilled. We can split the integral in two parts, resulting in 



f d''- l ki rdk . ( 



. .. ... , . . ... A A 2 d , //; 1 Am 2 

111= / ,.,./, / ^ln|*6 + — +#- .,. .^ ---^Afc 2 '— 



2 ' 2(d-l) V 4(d-l) 2 d-l 1 4 2{d-\) 



d c '- l kj fdk , (, 2 , m 2 , , 2 A / A 2 - , m 4 Am 2 



In + — +fcf- — rr -W -r; 1N9 - -3— rAff + ~T - ^73 ,J ■ (B3) 



(2;t) fl! - l 7 27t"'r n ' ' 2 ' 2(d-l) y4(rf-l) 2 d-l""* ' 4 2(d-l) 
Next, we can perform the integration over ko. In general, we can write, 

f dC ' k . r,2 \ f ,,2 ,2 \ f dd '% fd^o, ,2^ 



(B4) 
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as the second part is zero in dimensional regularization, (|99l l. Evaluating the integral over Icq gives, 



• jc — v / 



(27t)' / v y 7 (2t 
The second part is again zero in dimensional regularization so we obtain the following general formula, 

d d - l k: 



J {2%y K 



which we can apply to expression (IB3b 



III 



r d d -H 
J (2n) d - 



w-i 



A 2 



Afc? ■ 



Am 2 



\ 2 2{d-l) ]j 4(d-l) 2 d-1 > 4 2(d-l) 



id—1 



A 2 



{2n) d ^ \ 2 ""' 2(</-l) 
The next step will be to simplify this expression. If we define: 



4{d-l) 2 d-1 



d n m 4 
Akf- 



Am 2 



4 2(d-l) 



(B5) 
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(B7) 



2(d-l) ' 



A 2 



Am 2 



^-Afc 2 + — 

4(d-l) 2 d-l 1 4 2(d-l) 



we need to simplify the following expression, 



y/a + b+ ya — b . 



(B8) 



(B9) 



Notice that due to the constraint (1101b . a > for all values of kj. Let us first assume A > 0. For sufficiently small kf, b will be a 
positive real number, smaller than a. We may write 



((a + b) 1/2 + (a - fc) 1/2 ) 2 = ((a + fc) 1/2 ) 2 + ((a - fe) '/ 2 ) 2 + 2(a + b) 1/2 (a -b) x > 2 

= 2a + 2( fl 2 -/7 2 ) 1 / 2 , 

and taking the square of this equation results in 



Va + b+ \ / a-b= V2Va + \Ja 2 -b 2 



(BIO) 



(Bll) 



For kf larger than a certain value, the argument of the square root defining b will flip sign and b will become purely imaginary, 
i.e. b — ib' with b' > 0. The derivations ( IB 10b and ( IB 1 U remain valid, keeping in mind that (a + z'fo') always lies in the first 
quadrant of the complex plane, and (a — ib 1 ) in the fourth quadrant. 

For A < 0, we necessarily have that a > b > 0, which can be easily checked using the constraints ( I10U and (11031 ). Also now, 
(IB 101 ) and (IB lib go through, and we conclude that, given the original conditions ( 1101b and ( 1103b . we can always employ the 
equality in dBHI ). 

Using this formula, we can rewrite the integral III in the following form 



With the help of the following adapted Schwinger trick 



a+iVb 



■— lim 



2(d- 1) 



Jo Jo s 



kfiJki+A+m 2 



ds _, 2 f 
— — e 4j 

3/2 



-tA-sB 



(B12) 
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we can rewrite the square root 



III = 



%lJ (27i) rf - 1 \ z ^-i,Jo Jo s 3 / 2 V ' 



4- Km r a* r f / j^-^^ ) . ( bi 4 ) 



Therefore, we can now evaluate the integral over k, yielding 



{2Ti) d - 1 ' (2n) d - 1 Jo 



i ft 2 Y d/2 r(l) 



2 d-i n (d-i)/2 y4. s 'j r(^ T i ) 

If we insert equation dB 15b into equation ( IB 14b , we obtain, 

i r(i) r r s (d-3)/2 zihi^-siA+m 2 ) 

IE = = 7^- lim / dff z / ds ; -, 

^+i)/2r(«) z ^_i y Jo (t+2 S ) d 

= m I n/ , r( ; ^ lim f * e -^-^ ) f^' / ' (AW)/ V ( ''' 3)/2 (l+/)' rf , (B16) 

2 d/2 n (rf+l)/2 r (^) y Q t (d-\)/2-z J Q \ ' ' 

where in the last step, we have performed the substitution s' =2s/t. In this expression, we can switch the integral and the limit 
as it will turn out that the integral will converge (within the constraints ( 1101b and (11031 l). yielding, 

m = r* e - t (4-2^)rfc> e -MAW)/2 s ,(d-m {l+s , ) -d^ (Bn) 

2W+i)/2r(iB) Jo t d / 2 + l Jo K ' 

We recognize a Kummer function of the second kind, 

U(a,b,z) = 4- T dx^ x x a -\\+x) h - a - x 
T{a) Jo 

K 



sin(7t£ 



\F\{a;b\z) ( \F\(a — b + \;2 — b;z) 



T{a-b + \)V{b) V J T{a)T{2-b) 



(B18) 



where is \F\(a;b;z) a confluent hypergeometric function of the first kind. Comparing this Kummer function with expression 
dB17| ) we can write, 



-e 



2 d/2 % (d-i)/2 sin% d=i j fd/2+i y r(rf)r(-^=i) 



- -(A + m 2 ) i-^ — ,1V, • (B19) 

l 2 v V r(^ T i )r(^) y 



For the final integration, we recognize again a hypergeometric function, 

1 r + °° 

2 F 1 {a,b;c;z) = ^ / dfe^- 1 iFi(a,c,tz) , 

r{b) Jo 

resulting in 



(B20) 
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III 



2 d/2 n (<i-l)/2 sinn ^ 



- rf / 2 V 2 / / m 2 A \'/2 
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where this expression is defined in d dimensions. The next step is to replace d — > 4 — e, and to rewrite this expression in a series 
in e. Only in the first term, we do encounter a pole 1/e originating from r(— j). In the second term, there is no such pole and, 
therefore, we can immediately set d = 4: 



HI*, 



zF\ \ d, ~; '■ 



\(A+m 2 ) 



r(|) 



2 d/2 n (rf-l)/2 sinJt ^ 



>?i 2 A 

2 2(rf-l) 
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(m 2 +A) 5 /2 



1 7 w 2 



157I 2 



a2 A 
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2'2'm 2 -4 
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d=4 



For the first term, we have to expand in a series of e. For the benefit of the reader, we shall do this expansion in a structured way. 
We can distinguish 3 different parts. A pre-factor, the hypergeometric function and r(— |). Firstly, after some algebra, we can 
write the expanded pre-factor as, 



r(f) 
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2 <//2 7t (rf-l)/2 sinJt 



d-l 



(:i~nh)Y dl2 nd)n- d -^) 



8(47t) 2 



In I 



m 2 -4 



•ln(167t) 



)]) 



(B23) 



where we are working in the MS-scheme (later we shall convert /j to the MS scheme). Let us mention that in the calculation of 
this expansion, we have encountered the digamma function of —3 /2, i.e. V|l(— 3/2), which can be reduced to 



¥(-3/2) 



-21n2-y+8/3 : 



(B24) 



with the help of the following relations, 



¥(z) 



IM5 



v|/(z+l)-- 

z 



= -21n2-y. 

Secondly, we have to expand the hypergeometric function into a series in e. We find 
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We checked the explicit result for the expansion dB26b with the MATHEMATICA package HypExp 1271 |28J]. Finally, the expan- 
sion r(— i) in terms of £ reads, 



r(-^) = i 

v 2' e 



1-1 
4 2 



-0(e). 
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Taking the three previous expansions in e together, we find, 
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Notice that we have switched to the MS scheme. In summary, equation ( IB2U becomes, 
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